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Abstract 

We study the global existence issue for a three-dimensional Approximate Deconvolution Model 
with a vertical filter. We consider this model in a bounded cylindrical domain where we construct 
a unique global weak solution. The proof is based on a refinement of the energy method given 
by Berselli in [3]. 
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1 Introduction and notation 

The Large eddy simulation model (LES) with anisotropic regularization is introduced in [3] in 
order to study flows in bounded domains. 
Instead of choosing the classical filter 

A := / - {a\dl + a\d\ + a\df) , (1.1) 

the author in [3], considers a horizontal filter 

Ah := I - (ajdf + apl) . (1.2) 

This filter is less memory consuming than the classical one El [8]. Moreover, there is no need 
to introduce artificial boundary conditions for Helmholtz operator. However, the smoothing 
created by both of these filters is unnecessary strong. For this reason one can replace the 
classical Laplace operator by the Laplace operator with fractional regularization and seek for 
the limiting case where we can prove global existence and uniqueness of regular solutions. In 
a series of papers |10t HI [H [2], it was shown that the LES models which are derived by using 
instead of the operator (jl.ip . the operator Aq = I + a 2e (—A) e for suitable (large enough) 9 < 1, 
are still well-posed. 

In this paper, we consider LES models with fractional filter acting only in one variable 

A 3 , e :=I + af(-d 3 ) 2d , < 6 < 1. (1.3) 

In particular, we study the global existence and uniqueness of solutions to the vertical LES 
model on a bounded product domain of the type D = f2 x (— vr, 7r), where f2 is a smooth domain, 
with homogeneous Dirichlet boundary conditions on the lateral boundary d£l x (— 7r,7r), and 
with periodic boundary conditions in the vertical variable. For simplicity, we consider the 
domain D = {a;€ M 3 , x\ + x\ < d, —it < x 3 < 7r} with 2ir periodicity with respect to X3. More 
precisely, we consider the system of equations 
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dtw + V • w ® w — uAw + Vq = /, in D, t > 0, 
V • io = 0, in D, i > 0, 

Hanx(-7r,7r) =0, i > 0, 

W(XI,X2,X S ) = w(xi,X 2 , X 3 + 271"), t > 0, 

to (a;, 0) = too (a;) = W- 



(1.4) 



We recall that, in the above equations, the symbol " " denotes the vertical filter (|1.3p . applied 
component-by-component to the various tensor fields. Given that the filter is acting only in 
the vertical variable, it is possible to require the periodicity only in x 3 . Moreover, we consider 
the filtered function with homogeneous Dirichlet boundary conditions on the lateral boundary 
d£l x (—it, it). These boundary conditions of the filtered function are supposed to be the same 
as the unfiltered ones, in order to prevent from introducing artificial boundary conditions. 

The use of the vertical filter can be extended to a family of Approximate Deconvolution 
model. The deconvolution family including (|1.4p as the zeroth order case, is given by 



' d t w + V • D Ni o(w) <g> D N>e (w) - uAw + Vq = f, in D, t> 0, 
V • w = 0, in D, t>0, 
w\dnx(-TT,w) =0, t > 0, 
w(xi,X2,x 3 ) = w(xi,x 2 ,xs + 2ir), t > 0, 
, w(x,0) = w (x) = vfi, 

where -Djv.e is a deconvolution operator of order N and is given by 

N 



(1.5) 



D Nt0 = J2(i-&3})\ o<e<i. 



(1.6) 



i=0 



In order to state our main result pertaining to System (|1.5p . let us introduce our notation. 
We set: 

L 2 (D) := lv G L 2 (D) 3 , 2vr-periodic in x 3 , and j vdx 3 = o| , (1.7) 

endowed with the norm ||v||2 := / v ■ v dx which is the usual norm in L 2 (D) 3 . 

Jd 

Then, we define the following spaces that are commonly used in the study of the NSE: 



H := {v £ L 2 (D), such that V • v = and v ■ n = on dD) , 
V := {v £ H, such that Vv £ L 2 (D) and v = on dD} . 

Our main result is the following. 



(1.8) 
(1.9) 



Theorem 1.1. Assume f £ L 2 (0, T;H) be a divergence free function and vq £ H. let < 

N < oo be a given and fixed number and let 6 > i. Then problem M.5\) . including problem (1.4 ) 
when N = 0, has a unique regular weak solution. 

This result is one of the few results that consider the LES model on a bounded domain. It 
holds also true on the whole space M 3 and on the torus T3. Other a models, with partial filter, 
will be reported in a forthcoming paper. 

The remaining part of the paper is organized as follows. In the second section, we introduce 
several notations and spaces. We also prove auxiliary results which will be used in the third 
section. In the third section, we prove Theorem 11.11 Finally, a functional inequality is given in 
the appendix. 
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2 Preliminaries and auxiliary results 

Let 1 < p < +00 and 1 < q < +00. We denote by L q L p (D) = L q ((-vr, +vr); L p (n)) the space 

of functions g such that / I / \g{x\, X2, X2,)\ p dx\dx2) q ^ p dx^ ] < +00. In what follows, we 

denote by V/i the gradient operator with respect to the variables x\ and X2 ■ Of course, we have 
||Vfti?||2 < ||Vt>||2, for any »eV. 

Next, we will give some preliminary results which will play a important role in the proof of the 
main theorem. 

Lemma 2.1. For any u £ L 2 (D) 3 satisfying homogeneous Dirichlet boundary conditions on 
the boundary d£l x (— tt,tt) , such that V^u £ L 2 (L>) 3 , we have the following estimate 



II u\\ L 2 L 4 {D)3 < C\\u\\l \\V h u\\l (2.1) 

Proof. See in fTlf . 

Lemma 2.2. Let g be a 2ir-periodic and mean zero function with g{x) € H s ([—tt,tt]), and let 
s > ^, then we have the following inequality 

Iblk- <cNI^||<?||f s . (2.2) 

Proof. The proof is given in Appendix. 

Lemma 2.3. Let s > ^, and suppose that u and d s u £ L 2 (D) 3 , then there exists a constant 
C > such that 

i_J_ J_ 
IMIl-l2 (d)3 <C||u|| 2 2s ||d|u|| 2 2s . 



Proof. We first apply inequality (|2.2p in the vertical variable and then use the Holder 
inequality in the horizontal variable. We get, 



1 



sup I / |w(a;)| dx\dx2 I < I / sup |i*0e)| dx\dx2 



~ C (jn if ^ X ^ dx3 ) 2S (f \ d l u ( x )\ 2dx ^j ^ d Xl dx2^ (2 ' 3) 

< C\\Vu\\l~^ \\dlVu\\f . 
The following lemma will be useful for the estimate of the nonlinear term. 

Lemma 2.4. There exists a positive constant C > such that, for any s > \ and for any 
smooth enough divergence-free vector fields u, v, and w, the following estimates hold, 

00 • v)v,w)\ < c||u|||||Vu|||||v^||^^||a|v^|||^i|w|||||Vw|||, 

and 1111 1 

(ii) I ((u • V)v,w)\ < CllitlllllVnlllll^lllllVvlllUVwII^^II^VwIII 1 . 
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Proof, (i) By using lemma I2TT1 and lemma [23] we get 
\((u-V)v,w)\ 

- \\ u \\lIL\{D) 3 II ^ v \\l^LI(D) s \\ w \\lIL\{D) 3 



,1- 1 



(ii) We have 



< C||u|| 2 2 j|Vit|| 2 2 ||V«|| 2 23 ||<9f Vi>|| 2 2s ||w|| 2 2 1| Vif HI 



v ® u : Vw dx 

D 



\{(u-V)v,w)\ < 

thus, in order to prove the second inequality it is sufficient to swap the roles of v and w. 

2.1 The vertical filter and the vertical deconvolution operator 

Let v be a smooth function of the form v = X^fc 3 GZ\{o} Ck s ( Xl » x "i)^ kiX3 ■ The action of the verti- 
cal filter on v (x) = ^ c^ z (xi, X2)e lk ' AX ' A can be written as A^^^v) = Ae{k^)ck s {xi, X2)e lk3X3 , 

fc 3 ez\{o} fc 3 ez\{0} 
where the Fourier transform with respect to X3 of the vertical filter is given by 

A e (h) = (l + a 2e \h\ 2e y (2.4) 

Therefore, by using the Parseval's identity with respect to X3 we get, 



1 

k 3,9 v \\2 



2 v\\ 2 2 = \\v\\l + a 2e \\divf 2 = (A 3 , e v,v). (2.5) 



Next, we recall some properties of the vertical filter. These properties are proved in [3] for 
the horizontal filter and hold true for the vertical filter. 

Lemma 2.5. Let f be a smooth function and w a smooth, space-periodic with respect to x%, and 
divergence-free vector field defined on the domain D such that w ■ n = on dD. Let the filtering 
be defined by il.3\) . Then, the following identities hold true: 

(f,w) = (f,w), (2.6) 
w - a 2d dfw = w - a 2e dfw, (2.7) 

and 

(v -(w® w),A 3j9 u?j = (V • (w <g> to), w) = 0. (2.8) 

Proof. See in [3]. 

The deconvolution operator D^,e = Yld=oi^ ~ * s constructed by using the vertical 

filter with fractional regularization (|1.3p . For a fixed N > and for 9 = 1 , we recover a vertical 
operator form from the Van Cittert deconvolution operator (see [2] and [5] ) . A straightforward 
calculation yields 



Dn,9 [ Yl ^(xi,x 2 )e ik ^ 

K k 3 &Z\{0} 



26>lr, 1 20 \ N + 1S 



£ (W%| 2 *) U-( I fJ^ ) )c k3 ( Xl , X 2)e^ 



fc 3 GZ\{0} 
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Thus 

D N ,el 52 c fc8 (xi ) a:2)e ifc8Ba ) = £ V N>e (k 3 )c k3 ( Xl , x 2 )e ik ^, (2.10) 
\k 3 &\{0} J kez\{o} 

where we have for \ {0}, and 6 > 0, 

VoAh) = 1, (2.11) 

1 < V Nfi (k 3 ) < N + 1 for each JV > 0, (2.12) 

and T>Nfi(kz) < -4,3,0 for a fixed a > 0. (2-13) 

From the previous hypothesis, one can prove the following Lemma by adapting the results 
summarized in the isotropic case in [3]: 

Lemma 2.6. For 9 > 0, k 3 € Z \ {0} and /or eac/i ./V > 0, there exists a constant C > suc/i 
i/iai for all v sufficiently smooth we have 

\\v\\ 2 < \\D Nt g(v)\\ 2 < (JV + l)||w|| 2) (2.14) 

H| 2 < C\\D N>e (v) || 2 < C||A 3 ,^£)^( 1 ,)|| 2| (2.15) 

||A3,^D^(fj)|| 2 < H| 2 , (2-16) 

H£ + a 29 \\d e 3 v\\ 2 < ||Ag,^£>jt, 9 («)ll§. (2.17) 

3 Existence and uniqueness results 

In this section, we give a definition of what is called a regular weak solution to problem (|l,5p . 
Then, we give the proof of theorem 11.11 

Definition 3.1. Let f £ L 2 (0,T; H) be a divergence free function and Vq € H. For any 

< 6 < 1 and < N < oo, w is called a "regular" weak solution to problem il.5\) if the 
following properties are satisfied: 

w e C(0,T;H) r\L 2 (0,T;V), (3.1) 

d 9 3 w e C(0,T]H) nL 2 (0,T;V), (3.2) 

dfW G L 2 (0, T; V*), (3.3) 

and the velocity w fulfills 

fT 



j (d t w,ip) + (V- (D Nf g(w)®D Nt o(w)),(p) +v(Vw,Vcp) dt 

= [ T (J, tp) dt for all <p € L 2 (0, T; V). 
Jo 



(3.4) 



Moreover, 

w(0) = w . (3.5) 
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3.1 Proof of Theorem HTTP 



The proof of Theorem 11.11 follows the classical scheme. We start by constructing approximated 
solution w n via Galerkin method. Then, we seek for a priori estimates that are uniform with 
respect to n. Next, we take the limit in the equations after having used compactness properties. 
Finally, we show that the constructed solution is unique thanks to Gronwall's lemma [15j . 



Step 1 (Galerkin approximation). We denote by {f k3 }~^ = i the eigenfunctions of the Stokes 
operator on D, with Dirichlet boundary conditions on dD and with periodicity only with respect 
to X3. The explicit expression of these eigenfunctions can be found in [12]. These eigenfunctions 
are linear combinations of Wfc 3 (xi, X2)e lk3X3 , where k 3 S Z \ {0}, for certain families of smooth 
functions Wk 3 : O — >• R vanishing at dCl. 

We set 



w n (t,x)= 4 3 (^2,i)^ 3 ( 

fc 3 ez\{o},|fc 3 |<n 



X >- (3.6) 



Since we are dealing with real functions, one has to take suitable conjugation and alternatively 
use linear combinations of sines and cosines in the variable X3. We look for w n (t,x) which is 
determined by the following system of equations 

(d t w n , <p k3 ) + (^1%KP^KI, <^ 3 ) + v(Vw n , V</ 3 ) 

= (f,V k3 ), \k 3 \ = l,2,...,n, ( 3 - 7 ) 
w n (0) = P n (w ) = P n (vo), 

where P n denotes the projection operator over H n := Span(y> 1 , <p n ). Therefore, the classical 
Caratheodory theory [16] implies the short-time existence of solutions to (|3.7p . Next, we derive 
estimates on c n that is uniform w.r.t. n. These estimates imply that the solution of (|3,7p . 
constructed on a short time interval [0, T n [, exists for all t € [0, T]. 

Step 2 (A priori estimates) We need to derive an energy inequality for w n . This can be 
obtained by using A3 t gD]y t g(w n ) as a test function in (|3.7p . Thanks to the explicit expression 
of the eigenfunctions, and the properties of A3 and D^g, the quantity A^ ! gDj^ j Q(w n ) is a 
legitimate test function since it still belongs to H n . Standard manipulations and the use of 
Lemma 12.51 combined with the following identities 



(cW\ Aa.^Cto")) = \j t (ll4, e (^ n )ll2 + a 29 \\d e 3 Dl (w n )\\^j 
(-Aw n ,A 3 , e D Nte (w n )) = ( \\VDijw n )f 2 + WdiVD^Jw 



(3.8) 



.n\ ||2 
12 



|VA| e 4, e K l )ll2, 



(3.9) 



and 



(f,A 3 , e D N>e (w n )) = ^Dl e (f),Dl e (w n ^j , (3.10) 
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lead to the a priori estimate 



^114*^(^)11! + * /„ \\^ A !,9 D N,e( wn M >I« 



a 





D hW> D h( wn )) ^ + ^4^^0)111- 



(3.11) 



By using the Cauchy-Schwartz inequality, the Poincare inequality combined with the Young 
inequality and inequality (|2. 16|) . we conclude from ()3.11j) the following inequality 



(3.12) 



sup \\Al e D* Nd (w n )M + v \\VAI e D* Ne (w n Md S 

t€[0,T n [ JO 



< V 



2 , C(N + 1) 



112 



+ 



11/112 ds 



which immediately implies that the existence time is independent of n and it is possible to take 

J 1 — ipn 

We deduce from (f37T2l and (f2TT71) that 



sup 

te[o,T] 



>.W 



n\\2 



+ vJ*(\\Vw n \\ 2 2 + \\d!Vw n \\l) ds 



< v 



2 , C(N + 1) 



0112 



+ 



(3.13) 



" jo 

It follows from the above inequality that 

w n G L°°(0, T; H) n L 2 (0, T; V), 
dlw n G L°°(0, T; H) n L 2 (0, T; V). 

Thus, in one hand we get, 

Aw n G L 2 (0,T; V*), 
d e 3 Aw n G L 2 (0,T; V*). 

In the other hand, it follows from (|2.14p that 

D Ny6 (w n ) G L°°(0, T; H) n L 2 (0, T; V), 
^^(«> n ) G i°°(0, T; ff) n L 2 (0, T; V). 
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(3.14) 
(3.15) 

(3.16) 
(3.17) 



(3.18) 
(3.19) 



The inequality (|3.18|) allows us to find estimates on the nonlinear term V • {D^ ! e{w n ) <g> D^^(w n )). 
For that we use lemma l2~4l For all <p G V we have 



V • {D Nj6 (w n ) <g> Dj\[ t g(w n )), cpj < |(V • (r>jv,e(«? n ) <8) Djv,e(«> n )),^)| 
^CII^MIIallV^MIIallV^ll^llV^llf . 



(3.20) 



We also have 



e V-(D N , e {w n )®D N ,e( wn ))><P 



< 



V • (D Nj6 (w n ) ®D N>e (w n )),Aj tp 



(3.21) 
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Since ||V<9 3 V||f < ||V<^||f , ||VA|^||J 29 < \\V<p\\l 29 and || Vdf Af,^llf < l|V^||f we get 
that 

(v-(D Nt g(w n )®D Nt g(w n )),<p} \ <C\\D Ni g(w n )\\ 2 \\VD N> g(w n )\\ 2 \\V<ph, (3.22) 



and 



A| i0 V • (D Nfi (w n ) ® D N ,g(w n )),y 
Thus we obtain 



< C\\D N!0 (w n )\\ 2 \\VD N) o(w n )\\ 2 \\V<p\\ 2 . (3.23) 



V-(%K)®%K)) £l 2 (0,T;n (3-24) 

A| )0 V • (D Nt9 (w n ) ® D N , (w n )) € L 2 (0,T; V*). (3.25) 
From eqs. ([321), (|3"7Hj) . (f3T25|) . we also obtain that 

d t w n e L 2 (0,T;V*), (3.26) 

and 

a 4 A| e w n € L 2 (0, T; V*). (3.27) 

Step 3 (Limit n — > oo) It follows from the estimates ()3.14[) - (j3.27|) and the Aubin-Lions 
compactness lemma (see [13] for example) that there exists a not relabeled subsequence of w n 
and w such that 



w n - 


-* w 


weakly* in L°°(0,T;H), 


(3.28) 




D Nfi {w) 


weakly* in L°°(0, T; JT), 


(3.29) 


w n - 


w 


weakly in L 2 (0,T;F), 


(3.30) 


d e 3 w n - 


^d e 3 w 


weakly in L 2 (0,T;F), 


(3.31) 


D N , e (w n ) - 


^ D Nfi {w) 


weakly in L 2 (0,T;F), 


(3.32) 


d e 3 D N>e (w n ) - 


- d e 3 D N>e (w) 


weakly in L 2 (0,T;F), 


(3.33) 


d t w n - 


+ d t w 


weakly in L 2 (0,T;F*), 


(3.34) 


d t Aj >e w n - 


i 

- d t Al e w 


weakly in L 2 (0,T;F*), 


(3.35) 


w n - 




strongly in L 2 (0,T;H), 


(3.36) 


D N , e (w n ) - 


4 D Nfi (w) 


strongly in L 2 {0,T;H), 


(3.37) 



From (ET32|) and (13371) . it follows that 



V • (D N:d (w n ) ® Av,e(™ n )) V • (^(to) ® Av,e(™)) strongly in £*((), T; L 1 ^) 3 ), 

(3.38) 



Finally, since the sequence V • (D Ng (w n ) ® D Ne (w n )) is bounded in L 2 (0,T; V*), it converges 
weakly, up to a subsequence, to some x m ^ 2 (0> T\ V*). The previous result and the uniqueness 
of the limit allow us to claim that % = V • (D^ ; g(w) (g> D^^{w)). Consequently, 

V- {D N<g (w n ) ® D Nfi (w n )) ->> V • (D Ni g(w) ® D N<g (w)) weakly in L 2 (0, T; V*). (3.39) 
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The above established convergences are clearly sufficient to take the limit in (|3.7|) and conclude 
that the velocity w satisfies (|3.4p . Moreover, from (|3.30p and (|3,34p . one can deduce by a 
classical interpolation argument [U] that 

w G C(0, T; H). (3.40) 

Furthermore, from the strong continuity of w with respect to the time and with values in H, 
we deduce that w(0) = wq. 

Step 4 (Uniqueness) Next, we will show the continuous dependence of the solutions on the 
initial data and in particular the uniqueness. 

Let 6 > g and let {w\,qi) and (w2,q2) be any two solutions of (|1.5j) on the interval [0,T], with 
initial values tui(0) and 102(0). Let us denote by 5w = W2 — w\. We subtract the equations 
for Wi from the equations for W2. Then, we obtain 



d t Sw - vASw + V • (D Nj6 (w 2 ) <g> D Nfi (w 2 )) - V • (%(«i) ® = 0, (3.41) 

and 5w = at the initial time. 
1 

Applying A| to (|3.4ip we obtain 



A£ 9 d t 5w - uA^gASw + A|_ e V • (D Nj6 (w 2 ) <S> D N)0 (w 2 )) 

1 (3.42) 

-A 3 2 )0 V(-Av ie (™i) = 

1 

One can take A3 D^^ {5w) € L 2 (0,T;V) as test function in (|3.42p . Let us mention that, 

dtAjgdw e L 2 (0,T;V*) and A$ g D N>g (6w) € L 2 (0,T;V). Thus, by using Lions-Magenes 
Lemma [9] we have 

(d t Al e 5w,Al d n N , (6w)) v * )V = ^\\Al 6 Di/8w)\\l 

Using lemma [231 and the divergence free condition we get: 

— II^^C^IIa + ^IIVAl^CJtoJIll < [(DN t e(Sw) ■ V)D]y t0 (w 2 ), D N (5w)) . (3.43) 

We estimate now the right-hand side by using lemma [231 as follows, 
I ((D Ni0 (5w) ■ V)D N . e (w 2 ), D N0 (5w)) \ 

i_i J_ (3-44) 

< \\D Ni0 (6w)\\2\\VD Ni0 (6w)\\ 2 \\VD Ni0 (w 2 )\\ 2 26 \\8^VD Nfi {w 2 )\\^ 

Hence, by using the Young inequality combined with lemma [231 we obtain that there exists a 
constant C > that depends on TV and v such that 



I ((D N<0 (5w)V)D Nfi (w 2 ),D N j ) (5w)) 



< C\\Sw\\l\\Vw 2 \\l » ||flfViD2||| + ~||V<Hl2 



(3.45) 



By using (|2.17p we have 
1 d 



L± (\\ 5w f 2 + a^H^toIll) + (||V<H|! + c? e \\V5w\ 



(3.46) 
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From d33U), we get 

j t (\\5w\\l + a 2B \\dl5w\\t) + v (||V<Hll + « 20 ||W^||l) ^ ^ 



<C||<y«;||l||Vioa|ll 'II^Vm^HJ- 

2—- - 

Since HVW2II2 6 11^3 Vu^HI ^ -^([Oj^Dj we conclude by using Gronwall's inequality the con- 
tinuous dependence of the solutions on the initial data in the L°°([0, T], H) norm. In particular, 
if Swq = then 8w = and the solutions are unique for all t £ [0, T]. In addition, since T > 
is arbitrary chosen, this solution may be uniquely extended for all time. 
This finishes the proof of Theorem 11.11 

Acknowledgement I would like to thank Taoufik Hmidi for his valuable discussions, help 
and advices. 

Appendix: Proof of the functional estimate (12.21) 

Let g be a 27r-periodic and mean zero function with g(x) £ H s ([—ir, it]), and let s > ^. If we 

write g(x) as the sum of its Fourier series g(x) = gk& lkx , then, we can estimate ||<7||l°° 

fcez\{o} 

by 

< E \9k\- 

kez\{o} 

We then break up the sum into low and high wave-number components, 

< E \9k\ + E 

0<|fc|<K |fc|>« 

We now use the Cauchy-Schwarz inequality on each part, 

nu~< ( E i^i 2 ) ( E M + (ei^ 

\0<|fc|<K / \0<|A;|<k / \|fe|>K 

Since 

E 1 < Ck and E \ k \~ 2S ^ Ck~ 2s+1 , 

0<|fc|<K |fc|>K 

the above inequality becomes 

\\g\\L°° < C ^\\g\\ L 2 + K~ s+ ^\\g\\H» 

To make both terms on the right-hand side the same, we choose 

1 

WgWk* 

K = f-. 

\\9U2 

This yields the following estimate 
which is ([23]) . 
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